Coulomb versus spin-orbit interaction in few-electron carbon-nanotube quantum dots 



ON 
O 

o 



Andrea Secchi 1,2 and Massimo Rontani 2 

1 Dipartimento di Fisica, Universita degli Studi di Modena e Reggio Emilia, Via Campi 213/A, J h 1100 Modena, Italy 
2 CNR-INFM National Research Center S3, Via Campi 213/A, 4 1100 Modena, Italy 

(Dated: May 25, 2009) 

Few-electron states in carbon-nanotube quantum dots are studied by means of the configuration- 
interaction method. The peculiar non-interacting feature of the tunneling spectrum for two electrons, 
recently measured by Kuemmeth et al. [Nature 452, 448 (2008)], is explained by the splitting of 
a low-lying isospin multiplet due to spin-orbit interaction. Nevertheless, the strongly-interacting 
ground state forms a "Wigner molecule" made of electrons localized in space. Signatures of the 
electron molecule may be seen in tunneling spectra by varying the tunable dot confinement potential. 
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After almost two decades of research, carbon nan- 
otubes (CNs) 1] still provide a venue for the investi- 
gation of fundamental properties of interacting electron 
systems, such as Luttinger-liquid [H and Wigner-crystal 
y| behavior, Mott state Q , Kondo effect [Bj] and Andreev 
transport [6( in CN quantum dots (QDs). With respect 
to semiconductor QDs 0, S], low-screening, ultra-clean 
CN QDs appear to be ideal candidates [H for the real- 
ization of long-sought "Wigner molecules" (WMs) [U of 
strongly correlated electrons. These classical geometrical 
configurations of electrons localized in space are insensi- 
tive to the spin state of the system 0, Q • On the other 
hand, recently Kuemmeth et al. Q showed that orbital 
and spin degrees of freedom are entangled by strong spin- 
orbit interaction in single- wall CN QDs, disproving the 
popular shell model based on the fourfold degeneracy of 
QD energy levels 0,[lH. Intriguin gly, the two-electron 
tunneling spectrum measured in [j| was well explained 
by a non-interacting model with spin-orbit coupling. So 
far, it is unclear if the WM state [8| may survive to the 
perturbation induced by spin-orbit interaction. 

In this Letter we show that non-interacting features of 
the tunneling spectrum, due to spin-orbit coupling, co- 
exist with the strongly interacting nature of few-electron 
states, as seen from configuration interaction (CI, also 
known as exact diagonalization) calculations. Electrons 
in realistic dots form one-dimensional WMs, which may 
already have been observed in experiments [H, [§]. We 
predict that molecular signatures appear in the excita- 
tion spectrum by varying the QD confinement potential. 

We focus on a QD embedded in a semiconducting CN 
whose length scale, £qb, is smaller that the CN length. 
Hence, £qb is the relevant single-particle (SP) length and 
the effects of the CN boundaries may be neglected. With 
respect to previous calculations [TJ E2> E 0] > we as- 
sume the QD to be defined by an external gate poten- 
tial, slowly varying on the lattice scale, which we model 
as a one-dimensional harmonic oscillator (HO) of fre- 
quency wo- The quadratic potential is the low-energy 
generic form for a soft confinement [HI], setting £qy> — 
(h/m*LUo) 1 / 2 , where m* — ft 2 /3i?7 is the effective mass, 



R is the CN radius, and 7 = 0.54 eV-nm is the graphene 
7T-band parameter. SP states ip nTz (r) = NF n {x)<f) Tz (r) 
are obtained by the envelope-function modulation F n (x) 
of bulk states 4> T Jy r ) at t ne two non-equivalent minima of 
the lowest conduction band. The isospin index t z = +1 
(—1) labels valley K (K 7 ), F n (x) is the wave function of 
the nth HO excited state, and Af is a normalization fac- 
tor. Here ^(r) = exp(-iyT z /3R)[il>x,A{r)+T z ilJx,B(r)], 
where ^x,A{f) and ipx,B(r) are the Bloch tight-binding 
states for sublattices A and B, respectively, at point X — 
K (X = K') in the reciprocal space for t z = +1 (— 1). 
The isospin t z = +1 (—1) points to the (anti)clockwise 
rotation along the circumference coordinate y, perpen- 
dicular to the tube axis x. The interacting Hamilto- 
nian is H = Hsp + Vfw + Vbw, which includes the SP 
term H SP = Y,nT,<r z £ nTza*CnT x v,CnT z <r z , as well as the 
two-body terms for forward (FW) and backward (BW) 
Coulomb scattering processes Vfw and Vbw, respectively 
[3 14, [lj|. Here c nTz<Jz destroys an electron occupy- 
ing the SP orbital ip nTz (r) with spin a z . The SP energy 
includes the dominant term of spin-orbit interac- 



tion due to the CN curvature [1, [13, ll8| as well as the 
contribution of an axial magnetic field B [ll 



= s 



HO 



7 

- A S o S T zCr z 



*)■ (!) 



Here = r y/3R + hujo(n + 1/2), Aso is the spin-orbit 
coupling term, fj,B is the Bohr magneton, g* is the giro- 
magnetic factor. The two-body terms Vfw and Vbw are 

where 



CSo 



of the type V = 1/2 £ 4, 
a = (n, t z ) and V a ^s is the matrix element of the Ohno 
potential V(r - r') = U (l + e 2 \r - r'| 2 [/ 2 /e 4 )- 1/2 , 
which interpolates the two limits of Coulomb-like long- 
range and Hubbard-like short-range interactions (e is the 
relative dielectric constant, and Uq = 15 eV) [3]. In the 
envelope function approach [Hi ] the underlying graphene 
physics is buried into the precise form of V^/j^'s, which 
depend on both Bloch states 4> T X r ) and envelopes F n {x). 
We evaluate V a 7 s by conside ring explicitly the tight- 
binding expansion of 4> Tz ( r ) ( CI - ]XM)- FW and BW terms 
correspond to direct [{t z ,t z ) — > (t z ,t z )] and exchange 
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FIG. 1: (color online) Left: /j,(N) vs B for 1 < N < 3. 
Right: Main configurations contributing to states A and B. 
Blue (green) levels are HO states with isospin r z = +1 (— 1), 
whereas up- (down-) arrows stand for electron spin a z = +1 
( — 1). The percentages are the weights of the Slater determi- 
nants in the CI expansion of the wave function. 




t' z ] isospin scattering processes, 



The few-body problem is solved by means of the CI 
method [1, [2(|. We diagonalize H, which is a matrix 
in the basis of the Slater determinants obtained by 
filling with N electrons in all possible ways the thirty 
lowest-energy SP orbitals ip a (r). We obtain energies and 
wave functions of the many-body ground- and excited- 
states |\r/( n )), written as linear combinations of |$j)'s, 
|iT/(«)^ = J^c^ in each sector of the Fock space 
labeled by N, the z-component of the total spin S z , and 
the total parity under spatial inversion x — ► —x. 

To compare with tunneling spectra [§], we compute 
the chemical potential n{N) for a given value of B, 
M(i) (iV) = E 0(i) (N)-E (N-l), where E (N) (E. t (N)) is 
the energy of the many-body ground (ith excited) state 
with N electrons. The predicted /i(iV)'s may be con- 
verted into the gate voltages at which electrons tunnel 
into the QD [§] . We infer from [§] the inputs of CI calcu- 
lation, i.e., huj Q = 8 meV, R = 3.6 nm, A so = 1.24- 1(T 3 , 
g* = 2.14. The unknown value of e is a fit parameter for a 
small-gap semiconductor which may be strongly affected 
by the leads. By choosing e = 3.5 we obtain the curves 
of Fig. [TJ n(N) vs B for 1 < N < 3, which compare well 
with those of Fig. 3(a) in [jj. The plot quantitatively re- 
produces the dependence of fJ,(N) on B, specifically the 
kink of /i(2) at B c « 0.125 T and /x(3) at « 0.15 T [2l|. 
As a check of consistency with the experiment, we remark 
that the CI spacing between /x(2) and at B — 0, 17.8 
meV, agrees within 6% with the value estimated in Q. 

The kink of /x(2) at B c is due to the crossing between 
different ground states of two electrons, labeled A (B) for 
B < B c (B > B c ) [cf. Fig. H]. To gain insight into their 
nature, let us briefly consider also excited-state contribu- 
tions to Hi(N), which are shown in Fig.^a) for N = 1, 2. 
Note that Fig.fJ^a) perfectly matches Fig. 3(c) of 0. The 
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FIG. 2: (color online) Hi(N) vs B for N = 1, 2, for (a) htu = 8 
meV and (b) hujo — 4 meV. Inset: main configurations in the 
CI expansion of the wave function of state C. 



four lowest-energy SP levels /^i(l) are explicitly shown. 
Blue (green) lines indicate states with t z = 1 (—1), whose 
orbital magnetic moment is (anti)parallel to the CN axis, 
decreasing (increasing) its energy with B. The fourfold 
degeneracy of the levels (r 2 = ±1, a z — ±1) at B = 
is split by spin-orbit interaction, which entangles orbital 
and spin degrees of freedom 0]. This induces an isospin 
transition (t z = — 1 — > t z = 1) for the first excited non- 
interacting level /ii(l) at B m 0.125 T. Remarkably, this 
field is exactly the same as the critical value B c at which 
the A — > B transition occurs [cf. the vertical dashed line 
in Fig. [DJa)]. Because of this feature, it was argued in 
that the A — ► B transition may be solely explained 
in terms of spin-orbit interaction. The state A was sup- 
posed to be a single Slater determinant with the electrons 
in the two lowest spin-orbitals, (t z — 1, <j z — —1) and 
(t z = —1, a z = 1), whereas B is obtained by moving the 
electron from (r z = — 1, <y z = 1) to (t z = 1, cr z = 1). In 
this picture correlation effects are absent. We next show 
that A and B are instead strongly interacting states. 

The right panel of Fig. Q] shows the Slater deter- 
minants with the largest weights in the CI expansion 
of \A). The blue (green) ladders of levels depict the 
HO states for t z = +1 (—1), whereas arrows repre- 
sent spins. The main configuration, whose weight is 
53%, is the Slater determinant proposed in for the 
ground state. However, CI calculation shows that there 
are other three relevant determinants where also the ex- 
cited states of the HO are populated (Fig. [1]). This is 
due to the correlated character of \A), the strongest the 
impact of Coulomb interaction, the largest the mixing 
of determinants. Similarly, state B shown in the right 
panel of Fig. [T] is correlated as well. Besides, state B 
may be obtained from A by replacing the levels with 
t z = — 1 with those with t z = 1 (Fig. [T|). Hence, A 
and B belong to a isospin multiplet, only differing in 
the projections T z of the isospin. Here the isospin vec- 
tor is defined by f = ^i^^^^^ar^c^^ 
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where the components of <r are the Pauli matrices [e.g., 
T z = (n+i-n_i)/2 with n Tz = Y, n <? z 4t^c„^J- This 
is confirmed by extrapolating the behavior of A and B 
to the limit Aso ~~ ► (not shown). In this limit, A and 
B and other states at higher energies collapse into a six- 
fold multiplet, wich includes three spin triplet (T = 0, 
S = 1, S z = ±1,0) plus three isospin triplet (S = 0, 
T = 1, T z = ±1,0) states. In fact, the total wave func- 
tions must be odd whereas their (unique) orbital part 
is even under particle exchange. Therefore, except for 
a tiny residual splitting of w 2 /xeV due to BW inter- 
action, the A-B energy separation depends on spin-orbit 
interaction only. 

From Eq. (p} and the inspection of CI wave functions 
of Fig. [TJ it is clear that the energy splitting between 
A and B at zero field, /ii(2) — n(2) = 2Aso7/-R, is the 
same as that between one-electron levels, ^i(l) — /i(l) 
[cf. Fig.^a)]. Besides, n(N) depends on B through the 
spin and orbital magnetic dipole moments [cf. Eq. ([Tjl], 
which are linear in AS Z = S z (N) -S z (N-l) and AT Z = 
T Z (N)—T Z (N—1), respectively. It is immediate to verify 
that AS Z and AT Z are the same for both /i(2) and 
for B < B c (AS Z = 1/2, AT Z = -1/2) and B > B c 
(AS Z = 1/2, AT Z = 1/2). This explains why the critical 
value of the field for both /i(2) and ^i(l) is the same 
despite the correlated nature of A and B ; in agreement 
with the key experimental observation of [9( . 

The non-interacting feature of fj,(2) discussed above is 
not universal and may be affected by electron correlation. 
This is the case as one changes e.g. the QD potential, 
which is controlled in the laboratory by a capacitively 
coupled gate. This in turn changes the ratio of Coulomb 
matrix elements V to the SP spacing Hujq. Figure [2](b) 
shows the analogous plot of Fig. [2ja) for half the value 
of the confinement energy, hujQ = 4 meV. The pattern 
of n(2) vs B has now changed with respect to Fig. [Ha), 
due to the crossing between A and a new state, labeled 
C in the plot, occurring close to the origin, at B ss 0.02 
T. This critical value depends on the splitting Aca be- 
tween C and A at zero field, which in turn is sensitive to 
Coulomb correlation. 

The crucial role of Coulomb interaction in the QD is 
fully appreciated by considering the excitation spectrum. 
Figure [3] shows the two-electron excitation spectrum vs 
the dielectric constant e, which affects the relevance of 
correlation effects. The lower bound on the horizontal 
axis (e = 1.5) mimics low screening, typical of large- 
gap semiconducting CNs, whereas the upper bound (e = 
5000) may be regarded as the non-interacting limit. In 
the latter, all levels bunch into the HO levels, uniformly 
spaced by Hloq (= 8 meV). In the lowest set of levels 
(red lines) the two electrons occupy the nodeless orbitals 
in one of the two valleys, whereas the next set (black 
lines) is obtained by promoting one electron into the first 
HO excited state. The states of the first (second) set 
have even (odd) parity. Within each set, a residual fine 
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FIG. 3: (color online) Two-electron excitation energies vs e, 
reckoned from the ground state. Red (black) lines correspond 
to even (odd) parity states. 



structure survives, entirely due to spin-orbit interaction. 
As indicated in Fig. [H the total number of levels in the 
first (second) set is 6 (16), given by the possible ways to 
arrange the two electrons either in the same or in different 
valleys compatibly with Pauli's exclusion principle. 

As e is reduced in Fig. [3j Coulomb interaction alters 
the energy spectrum. In fact, the sixteen odd levels be- 
longing to the second set separate into two multiplets. A 
first sixfold multiplet is insensitive to e, whereas a second 
tenfold multiplet experiences a sudden energy drop. The 
former multiplet is associated to the collective motion of 
the center of mass (Kohn mode), which is decoupled from 
the relative motion and hence unaffected by Coulomb in- 
teraction [§] . The sixfold degeneracy of this odd multi- 
plet, lifted only by spin-orbit interaction, is the same as 
that of the lowest even multiplet, since they differ only in 
the excitation of the center of mass coordinate. The sec- 
ond odd multiplet is instead sensitive to e, i.e., Coulomb 
interaction. In fact, Aca becomes a small fraction of Hujq 
at the experimental value of e = 3.5 (Aca/^o ~ 0.06), 
whereas it vanishes at e = 1.5. As discussed below, such 
vanishing points to the formation of a WM, the state 
where Coulomb correlation localizes electrons in space to 
minimize their electrostatic energy 0, Q . 

The insets of Fig. 2] show the correlation function g(x) 
vs x for states A and C. g(x) is the probability of finding 
the two electrons at relative distance x — x± — X2 [nor- 
malized as J Q g(x)dx = 1, with g(x) — g{— x)]. As e 
decreases from e = 25 (bottom right inset) to the exper- 
imental value e = 3.5 (top left inset), the ground-state 
probability distribution <?a(^) (red curve) develops a well 
defined peak at xq — 1.7, showing that electrons local- 
ize in space and freeze their mutual distance, xq com- 
pares well with the equilibrium value xo.d = 1-6 of two 
point-like classical particles in the HO trap interacting 
via the Coulomb potential e 2 /e\x\. Note that crystal- 
lization proceeds by removing probability weight from 
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FIG. 4: (color online) dcA vs e. The dashed line is a linear 
fit to data in the range 5 < e < 15 with correlation coefficient 
0.999. Insets: 3a(x) and gc(x) vs x, for e — 3.5 (top left) and 
25 (bottom right). The length unit is ^qd- 



-B-value of the A — > C transition reported in Fig. &b) is 
a measure of the vanishing of Aca- As Hu>q is reduced, 
this critical field approaches zero, implying that the spin- 
polarized phase C may be induced with no energy cost. 
The latter behavior has been observed for hole WMs 0] . 

We have shown that spin-orbit and strong Coulomb in- 
teraction coexist in CN QDs, leading to the formation of 
WMs at experimentally attainable regimes. This insight 
into the entangled orbital and spin degrees of freedom is 
relevant for the all-electrical Q and -optical [22| manip- 
ulation of electron spins in CN-based devices. 
Note added. After the submission of this work, Wunsch 
reported similar results for a square- well QD [231 ]. 

We thank F. Manghi, E. Molinari, E. Andrei, G. Steele 
for stimulating discussions. This work is supported by 
INFM-CINECA Supercomputing Project 2008-2009. 



<7A at the origin. As localization is fully accomplished, 
§a{x — 0) — > 0, which is compatible with both even and 
odd states like A and C, respectively. In fact, gA(x) (red 
curve) and gc{x) (black curve) tend to coincide as well 
as Aqa — > as e decreases (cf. Fig. . Similarly, as the 
overlap between the probability weights of localized elec- 
trons is suppressed, exchange interaction is negligible and 
different values of S z (S z — for A and S z = — 1 for C) 
are admissible. In this limit states A and C, which only 
differ for now irrelevant quantum numbers like parity and 
spin, represent the same "classical" configuration. 

We exploit the progressive overlap between gA(x) and 
gc( x ) as e is reduced to characterize the transition to 
the WM state. In Fig. Q] we plot vs e the functional dis- 
tance dcA between states A and C, defined as dcA = 
J Q dx \gA{x) — gc{x)\. The semilog plot allows to iden- 
tify three separate regions, where dcA scales differently 
with e (there are no sharp transitions in finite-size sys- 
tems). For e > 20, dcA slowly tends to the upper bound 
1, since the location of the maximum of gA(x) approaches 
the origin whereas gc(%) has a node there. The large 
probability of finding two particles close to each other 
(x w 0) highlights the absence of a correlation hole in 
the ground state (cf. the magenta dot and related plot). 
In the crossover region, 5 < e < 20, dcA oc log e as shown 
by the linear fit in Fig. 0] (dashed line). Here a significant 
correlation hole rapidly forms in A as e decreases. The 
WM corresponds to e < 5, where dcA < 0.1 and it slowly 
decreases with e, as gA{x) and gc{x) overlap almost per- 
fectly. Remarkably, the observed case @ of e = 3.5 (blue 
dot in Fig. 2]) occurs in this region. 

The squeezing of the QD confinement potential via an 
external gate is a handle to drive Wigner crystallization, 
since the effect of lowering Itwq [Fig. [Hb)] is similar to 
that of decreasing e (Figs. [3] and 0}. In fact, the critical 
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